Heat capacity and pairing transition in nuclei 
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A simple model based on the canonical-ensemble theory is outlined for hot nuclei. The properties 
of the model are discussed with respect to the Fermi gas model and the breaking of Cooper pairs. 
The model describes well the experimental level density of deformed nuclei in various mass regions. 
The origin of the so-called S-shape of the heat capacity curve CV(T) is discussed. 



o ■ 
o ■ 

<n : 

Or 

<■ 

o ■ 
m ■ 



> 

OS 
(N 
O 

o 



X 
i 

=5 



X 



I. INTRODUCTION 

Nuclear structure at low excitation energy depends 
critically on the presence of Cooper pairs. Thermal and 
rotational breaking of these J = nucleon pairs gives 
abrupt structural changes, such as increased level density 
and rotational-spin alignments. These critical phenom- 
ena were early addressed in several theoretical papers 

A very exciting feature is the gradual reduction of 
pair correlations as function of temperature. Recently, 
Schiller et al. [Q reported the experimental critical tem- 
perature T c for the pairing transition. The findings were 
based on using the canonical heat capacity as thermome- 
ter. An S-shaped heat capacity as function of tempera- 
ture was observed in the 161 > 162 Dy and 171 > 172 Yb isotopes. 
Around T c ~ 0.5 MeV, a local maximum in the heat ca- 
pacity signals the breaking of Cooper pairs and quench- 
ing of the pair correlations. This property has also been 
observed in the 166 > 167 Er nuclei ||. 

Similar fingerprints have been obtained in various cal- 
culations. Finite-temperature Hartree Fock Bogoliubov 
calculations |@] for 164 Er give almost identical S-shape 
as observed for 166 Er. In relativistic Hartree Fock-BCS 
calculations the proton and neutron pairing gaps are 
seen to vanish around T ~ 0.4 - 0.5 MeV for ™,™Er. 
Furthermore, in Shell Model Monte Carlo simulations 
(SMMC) the heat capacities for iron isotopes show 
a pairing transition around temperatures of 0.7 MeV. 

The thermal breaking of a Cooper pair results in a ten- 
fold increase in number of available energy levels. In this 
process particles are thermally scattered on available sin- 
gle particle states, giving rise to increased entropy. Re- 
cently pCj^H l, it was shown that each thermal particle 
carries an entropy of ~ 1.7, a feature which is valid for 
midshell nuclei with mass number A > 40. 

The present work aims to present a simple model for 
hot nuclei that includes the main features found exper- 
imentally. In Sec. II the model is described within the 
canonical-ensemble theory, and in Sec. Ill some relevant 



model properties are discussed. The results are compared 
with recent theoretical and experimental data on various 
deformed nuclei in Sec. IV. Concluding remarks are given 
in Sec. V. 



II. MODEL 

The total partition function Z is described in the 
canonical ensemble, where thermal particle excitations, 
rotations, and vibrations are treated adiabatically ac- 
cording to 



Z(T) f=a Z par Z 10 tZ v ib. 



(1) 



Thermodynamical quantities such as the entropy S, the 
average excitation energy (E) and the heat capacity Cy 
can then be calculated from the Helmholtz free energy: 



F(T) = -T In Z(T), 



by 
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(E(T)) 
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(3) 
(4) 
(5) 



where the Boltzmann constant is set to unity (kg = 1) 
and T is measured in units of MeV. 

The theoretical basis for the particle partition function 
Z par was earlier presented in Ref. JTT[ , however, for com- 
pleteness we make a short outline here. The particle exci- 
tations are described by spin 1 /2 fermions scattered into 
doubly degenerated single-particle levels with equal en- 
ergy spacing e. There are infinitely many single-particle 
levels with one set of levels for protons and one set for 
neutrons. The partition function for one spin-up fermion 
is given by 
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zl = 5>xpHe/T) = r 



1 



exp(-e/T) 



(0) 



When two identical spin-up fermions are placed into the 
same single-particle level scheme, the summing is re- 
stricted due to the Pauli principle by 



E 



oo 



exp(-ie/T) exp(-je/T) 
=i+i 

1 exp(-e/r) 

1 - exp(-2e/T) 1 - exp(-e/T) 
t exp(-e/r) 
1 1 -exp(-2e/T)' 



(7) 



Generally, for n identical spin-up fermions the partition 
function is 



■A 



yl exp(-(n-l)e/r) 
l-exp(-ne/T) ' 



(8) 



on. The lowest partition functions for odd-mass (^1,3,5) 
and even-even systems (-^2,4) become: 

Zi = zi, 

Z 2 = 2z 2 , 

Z 3 = z 3 + Z1Z2, 

Z 4 = 2z A + (z 2 ) 2 , 

Z 5 = z 5 + z 3 z 2 + ziz A . (13) 

The partition functions for the odd-odd system can be 
constructed in a similar manner. The two-particle par- 
tition function is (z±) 2 , giving an entropy of twice the 
single particle entropy of the odd-mass system. Higher 
numbers of particles yield the following partition func- 
tions for the odd-odd case: 



Z2 = (zif, 
Z4 = 2z 3 z\, 
Z 6 = 2z 5 zi + (z 3 ) 2 . 



(14) 



We now introduce spin up and down for the fermions, 
and evaluate the corresponding partition function z n . For 
one fermion with spin up or down we obtain a degenera- 



tion of two, since z[ = z 



yl. 



z\ = z\ + z\ — 2z\. 



(9) 



If two fermions occupy the same single-particle level 
scheme, the partition function becomes more compli- 
cated. With the spins of the two fermions antiparallcl 
(to = toi + to 2 = 0), no Pauli blocking is present and 



z 2 (m = 0) 



T 1 
z[z{ 



(10) 



For the contribution with parallel spin, the Pauli princi- 
ple gives 



z 2 (m = ±1) = z\ + z\ 2 = 2z\. 



(11) 



It is now straightforward to evaluate z n for n fermions, 
allowing spin up and down. The five lowest partition 
functions are: 



z\ 
Z3 
Z5 



2z{ 
2z\ 
2z\ 



= 2z 



+ 2z 2 z^ , 

+ 2z 3 T z 1 + (z 2 )' 



2z\ 



2z\z\ 



^■z\z\. 



(12) 



The partition function for nuclei with protons and neu- 
trons may now be constructed from z n . In, e.g., the three- 
particle case of an odd-proton system the irv 2 and the tt 3 
partition functions have to be included. For an even-even 
two-particle system a n or a v pair can be broken, and so 



The creation of particles costs energy. In order to break 
one Cooper pair, the energy 2A is necessary. The pairing 
gap parameter is given by the empirical formula [ 12| 



A = 12A" 1/2 MeV. 



(15) 



The total particle partition functions of even-even (ee), 
odd (odd) and odd-odd (00) nuclei can then be expressed 
as 



Zpar — 1 



Z 2 e~ 2A ' T + Z 4 e- iA / T - 



Z;i d = Z 1 + Z 3 e- 2A / T + Z 5 e- iA / T - 
Z 2 + Z 4 e- 2A / T + Z 6 e- 4A / T - 



7OO 

^par 



(16) 



The number of active particles n (i.e. particles not bound 
in Cooper pairs) can be evaluated by 



odd 



I) 



{2Z 2 e- 2A ' T + AZ 4 e- iA ' T + . . .)/Z^, 
{Z X + 3Z 3 e~ 2A / T + hZ^ A ' T + . . .)/ZfJ : 
= (2Z 2 + AZ A e- 2A ' T + 6Z e e- 4A / T + . . .)/Z°° r . 

The rotational partition function is given by 

7=0,2,4... 

where the rotational parameter 



(17) 



(18) 



(19) 



is expressed by the rigid moment of inertia # r i g = 
2/5MR 2 ~ 0.0137.4 5 / 3 ft 2 MeV _1 , where M and R are 
the nuclear mass and radius, respectively. The vibra- 
tional partition function includes zero and one-phonon 
states: 
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Z vih = W v exp(-vhu vih /T), (20) 

where the multiplicity for the zero-phonon state is Wq = 
1. For one-phonon states, we take into accounts three 
vibrations (e.g. (3, 7 and octupole vibrations), giving 
multiplicity W\ — 3. Higher order phonon states are ne- 
glected, and the phonons are assumed to carry the same 
energy quantum %uj V ih. In both partition functions we 
have omitted the spin degeneracy (21 +1), since we are 
dealing with levels and not states. 



III. MODEL PROPERTIES 

Figure |l| shows the Helmholtz free energy F, the av- 
erage excitation energy (E) 1 the entropy S and the heat 
capacity Cy as functions of temperature. The model pa- 
rameters are taken from 162 Dy (see Table I) as a typical 
set applicable for rare earth nuclei. In order to calculate 
the thermodynamical quantities up to T ~ 1 MeV, at 
least 10 broken nucleon pairs have to be incorporated. 
The free energy F and the average excitation energy (E) 
behave smoothly as functions of temperature. Around 
T ~ 0.65 MeV the nuclei are excited to energies compa- 
rable to their respective neutron binding energies. The 
even-even, odd 1 and odd-odd systems have different exci- 
tation energies at one and the same temperature, where 
the even-even system requires the highest (E) value. 

The entropy S and heat capacity Cy represent first 
and second derivatives of F, and are thus more sensible 
to thermal changes. For the lowest temperatures the en- 
tropy difference is ~ 2 between the three mass systems. 
However, the entropy curves coincide for T > 0.6 MeV, 
displaying almost identical behavior. It is interesting to 
test if our model reproduces a Fermi gas with the absence 
of pairing at these temperatures 2 . The level density pa- 
rameter a for a Fermi gas with a level spacing e is given 

by@ 

7T 2 7T 2 7T 2 

0= y(ffp+9n) ~ Y 9= 3e' ^ 

Here, the single-particle level-density parameters for pro- 
tons and neutrons (g p and g n ) are assumed to be approx- 
imately equal. Inserting e = 0.13 MeV for 162 Dy, we ob- 
tain a = 25 MeV -1 . This is in exact agreement with the 
slope (a = dS/2dT) of the entropy curves for T > 0.6 
MeV. 



The lower right panel of Fig. [j] shows the typical S- 
shape of the heat capacity. As this shape is an impor- 
tant fingerprint for pairing transitions in nuclei dH0,fl > 
we will in the following focus on its origin. 

The contribution to Cy from collective excitations is 
negligible, and has no influence on the S shape. This is 
shown in Fig. for the even-even system, where the com- 
ponents of particle, rotational, and vibrational degrees 
of freedom are displayed. The collective contribution to 
the entropy S is small and fairly constant with increas- 
ing temperature. However, one should note that <SV t is 
in fact the main component at the lowest temperatures 
with T < 0.3 MeV. 

For the total entropy and heat capacity, also the effect 
of the number of pairs is shown in Fig. g. The Cy curves 
are identical for 6 and 10 pairs up to T ~ 0.7 MeV, but 
from there on the 6 pair system is exhausted and not able 
to absorb energy at the same rate as when more Cooper 
pairs are present. The figure shows that the S shape can 
be evaluated rather accurately without taking very many 
pairs into account. On the other hand, if too few pairs 
are included, one may easily misinterpret the shape of 
the Cy curve. For the 162 Dy mass region, 6 Cooper pairs 
are sufficient to determine the S shape. 

Figure || shows the heat capacity (upper panels) and 
the corresponding number of unpaired particles (lower 
panels) for 10 pairs of nucleons. In the left panels a re- 
alistic pairing gap parameter of A = 0.94 MeV has been 
chosen. We identify three temperature regions of inter- 
est: (i) The T ~ 0.2 - 0.8 MeV region with different S 
shapes for even-even, odd and odd-odd-systems, (ii) the 
Fermi gas regime with T ~ 0.8 — 1.0 MeV, exhibiting the 
same linear heat capacity for the three systems and (iii) 
the dramatical change for T > 1 MeV due to the final 
number of particles in the systems. 

The pronounced Cy maximum at T ^ 1.0—1.2 MeV is 
due the low number of available Cooper pairs; the lower 
left panel of Fig. || shows that ~ 15 nucleons are al- 
ready unpaired at this temperature. As T — *■ 00, we find 
Cy — > 20, 21 and 22 for the even-even, odd and odd-odd 
systems, respectively. This corresponds to the situation 
where all pairs are broken, giving a Cy-value equal to 
the total number of nucleons in the system 3 . Since a sys- 
tem with one particle in an infinite harmonic oscillator 
gives Cy — > 1 for T — > 00 and the effective Pauli blocking 
between the nucleons is negligible, we obtain Cy = n. 

In the right panels of Fig. ^, the pairing gap is reduced 
to A = e = 0.13 MeV. Now the S shape of the Cy curve 
can be seen to vanish in the A — > limit, telling that 



1 We use the abbreviation odd for odd-even and even-odd 
systems, since these systems are equivalent in our model. 

2 A simplified Fermi gas has S ~ 2aT + const. 

3 For the respective systems, we have 10 Cooper pairs plus 
0, 1 or 2 particles. 
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the S shape is connected with the pairing strength. Fur- 
thermore, we see that the heat capacity reaches a max- 
imum level of Cy ~ 35 at lower temperature T ~ 0.7 
MeV. At this point the remaining number of particles 
bound in Cooper pairs is low, and the depairing mecha- 
nism again looses its capability to create more heat ca- 
pacity. It might be surprising though that the number 
of depaired particles does not increase much faster with 
temperature. The answer to this is that in our model, 
not only the energy 2A is required to break up a pair, 
but also the unpaired particles have to be placed into 
the single-particle level scheme at some finite energy E, 
since the Pauli blocking prevents them to populate the 
lowest levels in the single-particle level scheme, which are 
all occupied by other unpaired nucleons. 

In Fig. |J an interpretation of the S shape is given for 
the even-even system. In the left panel, a "background" 
is subtracted from Cy, using a straight line through 0.3 
and 0.9 MeV of temperature. This line is intended to 
mimic the underlying heat capacity originating from a 
Fermi gas type of system. However, the gas properties 
arc strongly connected to the depairing process, and the 
subtracted peak in the lower left panel should be taken 
with care. 

The resulting peak shows some resemblance with the 
Schottky anomaly which describes a particle placed in a 
two level system |15J. However, the maximum heat ca- 
pacity of such a model is Cy ~ 0.45 at a temperature 
of 40% of the energy gap between the levels. For atomic 
nuclei, we may define a similar, but extended Schottky 
model: Either no pairs are broken at excitation energy 
E = 0, or one pair is broken at E = 2 A, or two pairs 
are broken at E = 4A and so on. This picture is an ex- 
treme simplification, and we include here only zero, one 
and two pairs. Thus, our Schottky-like partition function 
reads 



Zsch = 1 + W 2A e- 2A / T + W 4A e- 4A/T . 



(22) 



The multiplicities in front of the Boltzmann factors rep- 
resent the number of levels at E = 2 A and E = 4 A and 
are estimated from experimental data using W = p ■ SE, 
where p is the level density and SE is the energy window 
considered. 

In the right panels of Fig. [| Cy is evaluated by the use 
of Eq. (22) for 162 Dy. The p values are estimated from 
Table I of Ref. jn). With an energy window of 5E = 0.2 
MeV, we obtain Wsa = 22 and Waa = 556; numbers 
which are also consistent with the numbers of seniority 
S = 2 and 4 states in the model of Ref. Jl(| . The high 
multiplicities give a strength, width, and position of the 
peak that compare qualitatively well with the peak in 
the lower left panel. From this comparison, we interpret 



the local maximum of the S-shaped Cy curve at T ~ 0.5 
MeV as the point at which the depairing process is at 
the strongest. Hence, we define this point as the critical 
temperature T c for the pairing transition. 

To close this section, we comment on the weakening 
of the S shape when going from even-even to odd and 
odd-odd systems; an effect which is clearly seen in the 
lower right panel of Fig. [l]. The reason for the weakening 
can be understood from the entropy plot displayed in the 
right upper panel of Fig. |[ There, we saw that the en- 
tropy of each valence nucleon 4 is reduced from a value of 
S ~ 2 at T ~ 0.3 MeV to zero at T ~ 0.6 MeV. Since the 
corresponding heat capacity relates to TdS/dT, a nega- 
tive contribution to Cy appears and a quenching of the 
amplitude of the S-curve is apparent. 



IV. COMPARISON WITH DATA 

Our model is described within the canonical ensem- 
ble, while experimental data refer to the microcanonical 
ensemble. However, there are two ways to compare our 
model with experiments. With a known experimental 
level density p, the partition function can be constructed 
from 



Z(T) = J2SE lP (E l )e- E ^ T , 



(23) 



where Ei is the excitation energy and SEi are the energy 
bins used. From this partition function all thermody- 
namical quantities can be deduced, see e.g. Eqs. 
The drawback of this method, is that the level density 
function has to be known up to high excitation energies, 
typically E ~ 40 MeV. The other way is to evaluate the 
microcanonical level density from our canonical model. 
This can be performed by an inverse Laplace transforma- 
tion of Z. Using the saddle-point approximation (Fawler- 
Darwin method), we obtain p6[ 



P((E}) 



(24) 



The level density can then be compared with certain ex- 
perimental values, and complete knowledge on the level 
density up to high excitation energies is not necessary. 
One drawback using the saddle-point approximation is 
that the excitation energy (E) is a thermal average with 
a large standard variation of oe — T\fC ~ T\/2aT, with 
a being the level-density parameter. A second drawback 
is the approximation itself, which we have tested by a 
"forward-backward" Laplace transformation. In Fig. |5| 
we compare the original entropy and heat capacity, with 



Interpreted as the entropy gap between the various S-curves 
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the ones obtained by using Eq. ( |24| ) to obtain p, and then 
using Eq. (^3|) to obtain a new Z and its corresponding 
S and Cy. The comparison reveals a general smoothing 
for temperatures T < 0.5 — 0.6 MeV. The entropy is seen 
to be reproduced rather well, however, the heat capacity 
is more sensitive to the approximation. 

The model presented in Sec. II rests on the assumption 
that the single-particle level scheme can be approximated 
by equidistant levels. This assumption is never fulfilled 
in atomic nuclei, but within the canonical ensemble, the 
various deduced thermodynamical quantities are strongly 
smoothed with respect to excitation energy. Therefore, 
even with some non-uniformity, the model might still give 
realistic results. Probably, heavy and strongly deformed 
nuclei are the best candidates for our model. 

For test cases, we have chosen midshell nuclei around 
58 Fe, 106 Pd, 162 Dy and 234 U. All four mass regions re- 
veal a rather uniform Nilsson single-particle energy dis- 
tribution without large energy gaps. The pairing gap A 
and rotational parameter A ro t are calculated according 
to Eqs. ( |l5| ) and (p^). The vibrational energy quantum 
hujvib is taken as the energy at which the first vibrational 
states appear in the experimental level schemes of the re- 
spective mass region JlTj. The last parameter needed is 
the level-gap parameter e, which is expected to have a 
value in between the single-particle spacing e sp and the 
BCS quasiparticle spacing e qp = a/ (e sp — A) 2 + A 2 — A, 
where A is the Fermi level. 

In this work, e is chosen as a free parameter, de- 
termined from a fit to known experimental level densi- 
ties. For each nucleus we adopt two level density anchor 
points, as deduced in Ref. pfl . The lower anchor point 
is based on the counting of known discrete levels. This 
method is rather accurate, except for the odd-odd nuclei, 
where the number of levels might be several hundred per 
MeV, and thus difficult to measure. The other anchor 
point is based on average neutron-resonance spacing data 
at the neutron-binding energy. 

Figure^ shows anchor points and level densities 5 calcu- 
lated using Eq. (|24|). The parameters are listed in Table 
I, where e is adjusted to obtain the approximate slope of 
p in the log-plot. The agreement with the anchor points 
is good, and also the odd and odd-odd systems fall nicely 
into the systematics. For the 161 > 162 Dy isotopes, the ex- 
perimental level densities from Ref. Q] are shown as well. 
The calculations reveal good agreement with experiment 
as function of excitation energy. 

With the established parameters, we may now calcu- 
late the heat capacity Cy from our model and compare 
with experiments. Due to the imperfection of the saddle- 
point approximation, we take p((E)) from the saddle- 



point approximation (which is tuned to experiments in 
Fig. ||), and generate Cy from the corresponding Z func- 
tion of Eq. ( |23| ) . In the lower part of Fig. [7] the theoretical 
Cy curves are shown for 161 > 162 Dy. The experimental Cy 
curves, shown in the upper part, are based on experimen- 
tal level density data from Ref. || . Since the construction 
of Z requires data to much higher energies than experi- 
mentally known (E ~ 7 — 8 MeV) , the level density has 
to be extrapolated. Here, we have used the parametriza- 
tion of Egidy et al. p8[ , in accordance with our recent 
work B. The effective level density parameter (a ~ f7.5 
McV~^ has a major impact on the Cy curve for T > 0.5 
MeV. Figure ^ shows rather good agreement between ex- 
periment and model for T < 0.5 — 0.6 MeV. In particular 
the odd-even mass entropy difference is well reproduced. 
At higher temperatures the comparison is poor due to the 
arbitrary extrapolation of the experimental level densi- 
ties. 

In Fig. H the heat capacity Cy from our model, us- 
ing the parameters of Table I, is displayed for the four 
mass regions. The curves look similar; the differences are 
mainly the change in the scaling of the Cy and T axes 
for the various mass regions. The figure also includes the 
calculations on 58 - 59 Fe performed by Liu and Alhassad 
||. These data points show some discrepancies with the 
work of Rombouts et al. ||. However, both SMMC cal- 
culations obtain a critical temperature around 0.7 MeV, 
while we obtain a value around 1.2 MeV. Even when 
changing freely the number of particles, e and A in our 
model, we are not able to reproduce the SMMC results. 
In particular, the feature that Cy becomes constant 
or even drops || above T ~ 0.7 MeV is surprising. This 
effect might originate from shell gaps (which are not in- 
cluded in our model) or too few particles and/or orbitals 
considered in the SMMC calculations. It would be in- 
teresting to test if the SMMC calculations are capable 
of reproducing the anchor points for Fe in Fig. ^. Doing 
so, one should remember that our anchor points repre- 
sent densities of levels and not states, which are usually 
employed in SMMC calculations. 



V. CONCLUSIONS 

A simple model for hot nuclei has been outlined. The 
main properties of the model are determined by the 
pairing-gap parameter A and the energy-gap parameter 
e associated with an infinite single-particle level scheme 
for protons and neutrons. We have demonstrated how 
various thermodynamical quantities can be extracted. 

The model properties have been discussed with em- 



5 The experimental data are taken at a given excitation en- 
ergies E, while the calculations give p as function of average 
excitation energy (E). 
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phasize on the S shape of the heat capacity. This shape 
is strongly related to the bunch of newly created levels 
from the pair-breaking process at ~ 2A and ~ 4A of 
excitation energy. 

The model calculations are compared with experimen- 
tal data from the A ~ 58, 106, 162 and 234 mass regions. 
Using the saddle-point approximation with only one free 
parameter, the experimental level densities of even-even, 
odd and odd-odd systems are reproduced. 

The critical temperature of the pair transition for 
i6i,i62 Dy wag ca i cu i a ted to be T c ~ 0.5 MeV, in agree- 
ment with experiments. In the 58 - 59 Fe region we calcu- 
late T c ~ 1.2 MeV, which is significantly higher than the 
value of ~ 0.7 MeV obtained in Shell Model Monte Carlo 
simulations. Thus, further theoretical and experimental 
efforts are needed to understand the thermodynamics of 
these hot iron isotopes. 

We wish to acknowledge the support from the Norwe- 
gian Research Council (NFR). 
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Table I: Model parameters: A and A v - lg are taken from Eqs. ( |l5| ) and (JlJ), Tiuj v ^ from systematics, and e is tuned to 
fit the data of Fig. H. 



Nucleus 


A 


A r ig 


fi-Wvib 


e 




(MeV) 


(keV) 


(MeV) 


(MeV) 


58 Fe 


1.58 


42.0 


2.0 


0.80 


106 Pd 


1.17 


15.4 


1.4 


0.21 


162 Dy 


0.94 


7.6 


0.9 


0.13 


234 U 


0.78 


4.1 


0.8 


0.08 



6 




FIG. 1. Model calculations for nuclei around Dy. The four panels show the free energy F, the entropy S, the thermal 
excitation energy (E) and the heat capacity Cv as functions of temperature T. The same parameter set (Table I, 162 Dy) is 
used for even-even (solid lines), odd (dashed lines) and odd-odd systems (dashed-dotted lines). 
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FIG. 2. Entropy and heat capacity for Dy. The effects of the various degrees of freedom are displayed. The abrupt change 
when reducing the number of nucleon pairs from 10 to 6 pairs is evident. 
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FIG. 3. Heat capacity and number of unpaired nucleons for nuclei around Dy. In the right panels the pairing gap is 
reduced to A = e = 0.13 MeV. The calculations are shown for even-even (solid lines), odd (dashed lines) and odd-odd systems 
(dashed-dotted lines). 
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FIG. 4. Heat capacity for 162 Dy. The peaks in the lower panels are obtained by subtracting a linear background (dashed 
lines). In the right panels we have assumed a Schottky-like partition function, see text. 
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FIG. 5. Test on the saddle-point approximation. The solid lines are entropy and heat capacity obtained from our model for 
162 Dy. The dashed lines are calculated by making an inverse Laplace transformation with the saddle-point approximation and 
then back again with a Laplace transformation, see text. 
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FIG. 6. Calculated level densities of even-even (solid line), odd (dashed line) and odd-odd (dashed-dotted line) nuclei 
around 58 Fe, 106 Pd, 162 Dy and 234 U as function of average excitation energy (E). The open circles, squares and triangles 
are experimental level-density anchor points [11] extracted at certain excitation energies E. The solid circles for 161,162 Dy are 
experimental data points from Ref. [4] . 
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FIG. 7. Comparison between semiexperimental (upper panel) and theoretical (lower panel) heat capacity for 1 ' Dy. The 
semiexperimental values are taken from Ref. [4]. 
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FIG. 8. Calculated heat capacities for even-even (solid lines), odd (dashed lines) and odd-odd (dashed-dotted lines) nuclei 
around 58 Fe, 106 Pd, 162 Dy and 234 U. In addition, the heat capacities from SMMC simulations (scanned from Fig. 4 of Ref . [9]) 
are displayed as open triangles ( 58 Fe) and circles ( 59 Fe). 
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